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! Abstract. Maximum gravitational wave emission from binary stars in eccen- 

^ I trie orbits occurs near the periastron passage. We show that for a stationary 

>■ ■ distribution of binary neutron stars in the Galaxy, several high-eccentricity 



^ ■ systems with orbital periods in the range from tens of minutes to several 

, days should exist that emit broad gravitational- wave pulses in the frequency 

range 1-100 mHz. The space interferometer LISA could register the pulsed 
signal from these system at a signal-to-noise ratio level S/N > 5a/5 in the 
O ■ frequency range ~ 10~^ — 10~^ Hz during one-year observational time. Some 

^ ■ detection algorithms for such a signal are discussed. 

6 
u 

^ : 1 Introduction 



X 



The operation of the first broad-band laser interferometric antennas (LIGO, 
', VIRGO, GEO-600) is coming up The sensitivity of theses detectors in the 

frequency range 10-1000 Hz will be sufficient to observe astrophysical sources 
of gravitational waves (Braginsky 2000 and refs. therein). It is widely ac- 
cepted that compact binary stars consisting of neutrons stars (NS) and black 
holes (BH) are among the primary sources for these detectors (see Grishchuk 
et al. 2001 for a recent review). In the high-frequency band, which the first 
interferometers are designed for [|, the spiral-in phase of the binary coales- 
cences and their merging phase are to be observed, which can in principle 
be associated with spectacular astrophysical phenomena with an electromag- 
netic energy release of order of the NS binding energy (for example, with 
cosmic gamma-ray bursts, as was suggested by Blinnikov et al. 1984). 

^We remind that the sensitivity of the ground-based detectors is Hmited at low fre- 
quency by seismic noises 



1 



In the case of two material points of equal mass in a circular orbit with period 
P, the coalescence time on the quadrupole approximation is approximately 

to ~ lQ%P/[ls]f' [s] 

i.e. about 20 minutes for the initial P — 0.2 s corresponding to a GW 
frequency of 10 Hz at the low-frequency sensitivity boundary of the ground- 
based interferometers. The signal has a quasi-periodic character with grad- 
ually increasing amplitude and frequency as the stars approach each other 
(the so-called "chirp" signal). The waveform of this signal can be calculated 
with a very high accuracy so that a powerful optimal filtering technique (well 
known in radiophysics, e.g.. Tikhonov 1966) can be applied to extract the 
signal out from the instrumental noise (see Grishchuk at al. 2001 for more 
references). 

At lower frequencies, the searching for GWs is possible using space antennas. 
The most advance project to date is the LISA interferometer (Bender et al. 
2000). The planned sensitivity range for this detector lies within 10~^ — 10~^ 
Hz. Main sources for LISA include massive binary BH coalescences in remote 
galaxies, individual close binaries in our Galaxy, and, possibly, a relic GW 
background originated in the very early Universe (Grishchuk et al. 2001). 
Galactic close compact binaries (white dwarfs, NS, BH) mainly compose a 
stochastic GW noise at frequencies below 1 niHz. The binary stochastic GW 
background and the possibility of detecting individual binaries with compact 
objects by LISA was recently reconsidered by Nelemans et al. (2001). 

In this letter, we study GW-signals from high-eccentricity binary NS, which 
appear as a sequence of broad-band pulses with a period from tens of minutes 
to several hours in the frequency range 10~^ — 10~^ Hz. These pulses form 
during the periastron passages of binary NS, where most GW are emitted. 
For a stationary galactic binary NS coalescence rate of 10~^ events per year, 
the number of binary NS that can produce such pulses at frequencies 1-100 
niHz at the LISA signal-to-noise level 5\/5 is estimated to be around 7 for 
one-year observation time. 
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2 Stationary distribution function of binary neu- 
tron stars in the Galaxy 



Binary NS (and BH) are natural products of evolution of massive stars. In 
the case of NS, core-collapse type II or Ib/c supernova explosions occur twice 
during the evolution of the original massive binary. It is well known (Blaauw 
1961) that the sudden mass-loss from one component of a circular binary 
system due to supernova explosion makes the orbit non-circular with an 
eccentricity e = AM/M', where AM is the mass lost and M' is the total mass 
of the system after the explosion. If an additional kick velocity is imparted 
to the newborn NS (which largely follows from radio pulsar space velocity 
observations, e.g. Lyne and Lorimer 1994, Arzoumanian et al. 2001, or is 
directly indicated by the non-coaxiality between the NS spin and the orbital 
angular momentum in some binary pulsars, e.g. PSR B1259-63, Prokhorov 
and Postnov 1997, PSR J2051-0827, Doroshenko O. et. al 2001), the resulting 
binary eccentricity become even larger on average. 

In addition to having a sufficiently large amplitude, a GW signal should be 
frequent enough to have chances to be detected during a reasonable obser- 
vational time (usually taken to be one year). At present, a good theoretical 
estimate of the galactic binary NS coalescence rate is 10~^ — 10~^ per year 
(Lipunov et al. 1997, Yungelson and Portegies Zwart 1998 and later papers, 
see Grishchuk et al. 2001 for a more detailed discussion). The main theoret- 
ical uncertainty pertinent to such estimates is due to unknown anisotropy of 
the stellar core collapse. Since only close system that evolve rapidly enough 
can coalesce due to GW emission over the galactic age ~ 10^*^ years, the 
galactic birthrate of binary NS must be higher. So it is expected that under 
stationary conditions N — 10^ binary NS can simultaneously exist in the 
Galaxy. 

Assuming that GW emission is the only mechanism for the orbit decay (which 
is a very good approximation in sufficiently low-dense stellar fields; in dense 
stellar systems dynamical interactions should be taken into account) and 
that a stationary star formation has really taken place in the Galaxy over 
past 5 billion years, one can calculate a binary NS distribution function over 
orbital parameters (orbital large semi-axis a and eccentricity e), which is 
normalized as N — J dN{a,e)dade. 

This problem was solved for model initial parameters by Buitrago et al. 
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(1994). In our previous paper (Ignatiev et al. 2001), the distribution func- 
tion dN{a, e) was found for the initial orbital parameters of binary NS com- 
puted using the population synthesis code "Scenario Machine" (Lipunov et 
al. 1996). In the parameter region where the binary coalescence occurs over 
the time interval shorter than the galactic age, the distribution function be- 
comes stationary. A fragment of the stationary distribution calculated using 
the Maxwellian distribution of the kick velocity with a mean value of 200 
km/s is shown in Fig. 1. The normalization is made for the galactic binary 
NS coalescence rate 10~^. Clearly, the function sharply increases at low or- 
bital frequencies (large orbital periods). The shape of the initial distribution 
function has most effect exactly in this region. 



3 GW-signal from binary NS at periastron pas- 
sage 



A binary star with the components of mass Mi and M2 in a non-circular 
orbit emits gravitational waves in a wide frequency range on harmonics to the 
orbital Keplerian frequency (27ri/x)^ = G{Mi + M2)/a^ (Peters and Mathews 
1963). The main energy is emitted when the stars approach the periastron 
distance = a(l — e). For large eccentricities, the duration of the GW 
pulse is approximately r ~ ap/vp (vp is the velocity at the periastron) and 
the maximum of radiation is emitted in the harmonic whose number n ~ 
(tuk)'^ ~ (1 — e)~^/^ in accordance with the 3d Kepler's law. 

Let us calculate more precisely the time of the most energy emission during 
the periastron passage. We shall use the well known expression for GW 
emission from a binary in elliptical orbit averaged over the binary orbit 
orientation (Peters and Mathews 1963, Landau and Lifshitz 1988) 

-— = 15^5,5(1 _ e2)5 H I + e cos 0)^[12(1 + e cos 0)^ + sm^ 0] . (1) 

Here and below G is the Newton gravitation constant, c is the speed of 
light, (f) is the polar angle counted from the periastron Next, we define the 
periastron GW pulse duration r as the time of emission of 90% of the energy 
released over one orbital period P. Let us change integration over time t by 
integration over the angle 



;i_e2)3/2 V a 



l + ecos0)2 G{Mi + M< 
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(we assume that GW-induced change in the orbital energy and angular mo- 
mentum is small and does not affect the Keplerian parameters of the orbit; 
this is justified for orbital frequencies considered). The dependence of the 
polar angle (690 that determines the corresponding part of the orbit is shown 
in Fig. 2. The plot r/P versus eccentricity is also shown in Fig. 2. Starting 
from e ~ 0.5 this angle is about 60 degrees and slowly changes with eccen- 
tricity (which is intuitively clear: at large e the focal parameter of ellipse 
p and the periastron distance are both proportional to (1 — e) and the 
periastron parts of orbits with different e are similar). As also follows from 
Fig. 2, already at e ~ 0.3 90% of the energy is released over 1/3 of the 
orbital period. Below we shall take e = 0.3 as a low boundary of binary NS 
eccentricities which can produce broad GW pulses. 

Thus, a binary NS in an elliptical orbit emits a GW signal in a wide fre- 
quency range Az/^ ~ l/r. The pulses are characterized by the duration r, 
the recurrent period P, and by some amplitude h. The latter cane be es- 
timated as a maximal dimensionless GW amplitude h^ax at the periastron. 
Averaging the squares of the field amplitudes h+ and h^. at the periastron 
{(f) = 0) (Peters and Mathews 1963, Moreno-Garrido et al. 1995) over the 
binary star orbit orientation with respect to the line of sight, we arrive at 



V 5 (^a{l-e^)r^ 3 3 3 ^ ' 

where r is the distance to the system. In the limit e — > 1, this maximum 
amplitude exceeds by ^28/3 ~ 3 times the amplitude of GW in a circular 
orbit with radius equal to the focal parameter p, and by ^-^ 1.5 times that with 
the periastron distance. We also note that the mean amplitude y^Ti^r-i-"^ 
for the case of rgo changes with increasing eccentricity at the endpoints of 
the corresponding time intervals by no more that 2 times, and by no more 
than 4.5 times for the case of rgg (Fig. 2, the bottom panel). 

For numerical estimates we should know the distance to the source. Since 
during a binary NS formation the barycenter of the system can acquire an 
appreciable space velocity, the galactic subsystem of binary NS has a larger 
size than that of ordinary stellar components. Quantitative calculations were 
performed by different authors (e.g. Bulik et al. 1999), so we shall use their 
results. For example, for the kick velocity 200 km/s the mean distance to 
a binary NS from Sun is about 12.15 kpc. For numerical estimates we shall 
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assume all binary NS stars to lie at this average distance. 



4 Detecting broadband GW signals 



To calculate the expected detection rate of GW signals considered it is nec- 
essary to precise methods of their possible detection. First, let us choose 
the frequency range in which the signals are to be searched for. It is known 
(see Grishchuk et al 2001, Nelemans et al 2001) that the detection of GW 
signals in the LISA frequency band up to ~ 1 mHz is hampere(^ by the un- 
avoidable presence of a stochastic noise formed by unresolved galactic binary 
white dwarfs (the stochastic background associated with unresolved galactic 
binary NS even with account of non-circular systems is several times lower 
due to much lower amount of binary NS than white dwarfs, see Ignatiev et 
al. 2001). So in what follows we shall consider only the frequency range 
10-3 - 10-1 Hz. 



Binary NS in elliptical orbits whose GW emission at periastron falls within 
this frequency range have orbital periods P from tens of minutes to several 
days (in the case of very large eccentricities). That is, during one year ob- 
servation time (T ^ 3 X 10^ s) the signal from one system will consist of a 
periodic sequence with k = T/P broad pulses with an arbitrary (from the 
beginning of the observation) phase. The direction to the source is also un- 
known. Clearly, there should be more than one such systems (cf. Nelemans 
et al. 2001). So from the viewpoint of detection the total signal is a super- 
position of periodic pulse sequences Si(t) with arbitrary phases and different 
periods s{t) = si(t) + S2{t) + . . . + s„(t). We stress that in the frequency 
range considered such a signal does not form a continuous background with 
the detector frequency resolution Au = 3 x 10"^ Hz (Ignatiev et al. 2001) 
and the planned LISA sensitivity Sn- 

Consider now an ideal case where there is only one sequence of such pulses. 
The signal-to-noise ratio when observing one broad pulse with the duration 
r by a detector with a sensitivity Sn [Hz-^/^] is 

(S/N), = (3) 

where Auj, = 1/r is the pulse frequency width. For a periodic sequence of k 
^But is not totally excluded, see the resent paper by Hellings (2001) 
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identical pulses we obtain 

{S/N'k = {S/N),Vk = ^^./- (4) 




where Au = 1/T = 3 x 10 ^ Hz is the maximum frequency resolution over 
the continuous observation time. 



As an example, consider a pulsar similar to the Hulse-Taylor's one PSR 
B1913+16 with orbital parameters (Taylor and Weisberg 1989) a = 1.95 x 
10^^ cm, P = 27907 s, e = 0.617, the equal component masses Mi = 
M2 = IAMq, located ~ 5 kpc away. Averaging over the binary orienta- 
tion, we find hmax ~ 10~^^. With the eccentricity 0.617 the pulse duration is 
Tgo ~ 0.15P = 4186 s, Az/fc ^ 0.23 mHz, the frequency of the maximal (4th) 
harmonic is ~ 0.18 mHz, and the signal-to-noise ratio assessed using Eqn (Q) 
is only 0.5 (without taking into account the binary white dwarf noise in this 
range). Clearly, such sources are of no interest for us. 

For reliable detection of a pulse signal coming from unknown direction one 
usually requires S/N > where the factor y/b accounts for unknown co- 
ordinates of the source (Thorne 1987). With such a high detection threshold 
the integration of the stationary distribution function over the region where 
binary NS systems giving such a high amplitude reside is ~ 7. Clearly, the 
lower the detection threshold the larger the number of the system (for exam- 
ple, it is about 15 for S/N = 5). Increasing the pulse duration (e.g. taking 
Tgg instead of Tgo) could formally somewhat increase the number of detected 
systems due to decreasing the frequency band Auf, in Eqn (||), however in this 
case the systems should be chosen starting from a higher eccentricity (~ 0.5) 
in order that the pulse duration be less than 1/3 of the orbital period. So in 
fact the corresponding region of integration of the distribution function gets 
narrower and the number of detections decreases. 



The distribution of binary NS for the stationary distribution function under 
consideration over the pulse widths r and orbital periods P is shown in Fig. 
||. The maximum probability is to detect pulses about 100 s in duration 
from systems with orbital period near 1000 s. The mean eccentricity of such 
systems, as follows from Fig. 2, the middle panel, is about 0.5. 

Thus, actually several systems with sufficiently high amplitudes simultane- 
ously contribute to the total signal. Here we meet the typical problem of 
parameter estimation of a quasi-determined signal in the background noises. 
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For Gaussian noises and the a priori known number of sources this prob- 
lem can be solved exactly using the well known methods of signal detection 
in radiopliysics. The basic principles of constructing optimal detectors for 
quasi-determined signals with Gaussian noises are collected in the Appendix. 



5 Discussion 



1. Stationarity of binary NS distribution in the Galaxy over orbital parame- 
ters is one of the main assumptions used. This assumption is justified when 
the characteristic evolutionary time scale of the systems is much shorter than 
the time of significant variations of the initial source distribution function. 
We assume that the star formation rate in the Galaxy can be considered 
constant at least over the last 5 billion years (a different assumptions is used 
in the paper by Nelemans et al. (2001)). The evolutionary time scale of 
compact binaries in the frequency range considered 10~^ — 10~^ Hz is much 
shorter than this time interval so the stationarity can be established. The 
increase of the star formation rate in the past leads to some increase of the 
value of the distribution function, especially at low frequencies. However, 
parameters of binary NS formation (especially the value of the kick velocity 
and its distribution) has more effect on the actual binary NS distribution. 

2. Normalization of the distribution function to provide the binary NS co- 
alescence rate 10~^ per year is a good estimate which is close to the upper 
accepted boundary of such evaluations (Grishchuk et al. 2001 and refer- 
ences therein). Changing this normalization would lead to the corresponding 
change in the resulting numbers. 

3. We do not consider a more numerous potential class of sources like a NS + 
a massive white dwarf. Systems in circular orbits are outside our scope, and 
systems NS+massive WD in elliptical orbits (studied for example by Porte- 
gies Zwart and Yungelson (1999)) in the frequency range of interest are not 
distinguishable from double NS systems. A detailed analysis of formation 
of binary systems comprising an old white dwarf and a young NS, thereof 
representatives can be binary radio pulsars PSR B2303 46 and PSR .11141- 
6545, has recently been performed by Tauris and Sennels (2000). According 
to these authors, the (model-dependent) formation rate of such systems in 
the Galaxy can substantially exceed that of binary NS (see also Brown et al. 
2001). We also omit compact binary NS+BH and BH+BH systems due to 
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their small number in comparison with binary NS. Anyway, the possible con- 
tribution of systems like NS ^ WD, NS l-BH, BH+BH in non-circular orbits 
will increase the number of pulses analyzed. 



6 Conclusion 

We have analyzed the broad pulse GW signal produced by compact binary 
stars during the periastron passage. In the particular case of a stationary 
distribution of binary NS in the Galaxy over orbital periods and eccentricities 
normalized to the galactic coalescence rate 10"'' per year, in the frequency 
range 10~^ — 10~^ Hz, one can expect to detect with the LISA space in- 
terferometer ~ 7 individual sources at the signal-to-noise ratio level 5-\/5. 
Allowing for other types of compact binaries (NS with white dwarfs or black 
holes) can increase this number. The GW signal from each system represents 
a sequence of broad pulses with equal amplitude. The most plausible width 
of the pulses is around 100 s, the period (the orbital period of the underly- 
ing binary system) is around 1000 s. From th viewpoint of signal detection 
theory, some optimal detection algorithms of such a quasi-determined signal 
in the background of additive noise are discussed. 

The authors acknowledge L.R.Yungelson for constructive notes. The work 
is partially supported by the RFBR grants 00-02-17164, 00-02-17884a, and 
01-lo-88310m. KAP thanks the MPA fiir Astrophysik (Garching) for hospi- 
tality. 



7 Appendix 

Let 

x{t) ^ s{t,ei,...,eN) + n{t) , O^t^T 

be an additive mixture of a quasi-determined signalf5'(i) 

N 

S{t, Gi, Qn) = J2 Q^"^^ W ' O^i^N 

1=1 

where = (0, 1) is the detectability parameters, and a Gaussian noise n{t) 
with the known correlation function /^^(t). 



9 



A. Detectability parameters are known. The hypothesis that the signal is 
present is accepted if the following condition is met (e.g. Sosulin 1992) 

N N N 



2 

1=1 i=l k=l 



Here 



li = y" x{t)ui{t)dt , Pik^ J Si{t)uk{t)dt , 



Y ^ {¥,,... ,Y^) . 

where the fiducial signal Ui{t), i — 1, N is the solution of the type I Fredholm 
equation 

T 

J Knit - T)Uk{T)dT = Skit) , O^t^T . 


The threshold level depends upon the false alarm probability a (Neumann- 
Pierson test). 

The analysis of Eqn (Al) shows that for known detectability parameters 
the optimal detector represents an A^-channel linear system. The princi- 
pal element of an individual system is a correlation detector that from the 
independent variable Y^, i — 1,N (in the filtration variant the correlation 
detector is substituted by an optimal filter). 

B. Detectability parameters are unknown. Bayesian algorithm of detec- 
tion. Consider the detectability parameters ©j to be independent discrete 
random values with known distribution 



p{e, = i} = p,, p{e, = o} = i-p, , t = i,N. 

It can be shown that under the condition of parametric apriori uncertainty 
the hypothesis that the quasi-determined signal S{t,Qi, ...,Qn) is present 
can be accepted if 

Z{Y,e^B,...,ONB)>Co^ , {A2) 

where G^^ is a Bayesian estimate of the unknown parameter with a quadratic 
payment function 

^iB J ...J e,exp[z{Y,Q,,...,Q^)^Wpr{Qi,...,QN)de,...deN , 



N 



{A3) 
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N 
i=l 

is an apriori joint probability density of random variables ©i, ©jv- 

C. Detectability parameters are unknown. Non-Bayesian detection algo- 
rithm. For non-parametric apriori uncertainty the detectability parameters 
©j are considered as non-random. When constructing a detection device, 
parameters ©j are substituted by maximum likelihood estimates 0j, which 
leads to the following discrimination rule: 

z(y,©i,...,©jv) >c„ . (A4) 

The maximum likelihood estimates ©j are chosen such that 

Z(f,0i,...,0^) > max ©i, ©^^), (A5) 

The analysis of Eqn (A2-A5) indicates that: 

1) For detecting quasi-determined signal S{t, ©i, ©jv) the detector should 
contain an additional unit, the estimator of unknown detectability parame- 
ters ©1, ©JV- 

2) Detection algorithms (A2) and (A5) prove to be non-linear: 0j = 0j(F), © 

3) For a non-Bayesian detection algorithm, the block of estimates of unknown 
detectability parameters represents a multi-channel linear system (see (A5)). 
The output signal from this systems is used in the matching scheme. 

4) The detection of the quasi-determined signal ©i, ©jv) using algo- 
rithms (A2) and (A4) can be considered as a joint detection and measurement 
of parameters ©j of the original signal. 

5) The structure of the Bayesian estimation block is essentially non-linear so 
such system are not widely used in practice. 

The threshold level in Eqn (A2) and (A4) is determined by the following 
conditions 

P{Z{Y, ©IB, ©jvb) > Ca|0i=02=...=ei=o} = a , I 

... } m 

P{Z{Y, ©1, ©^) > Ca\e,=e,=...=e,=o} = « , I 



where in Eqn (A6) ©is and ©j are pseudo-estimates of unknown detectability 
parameters in the absence of the signal. 

Since algorithms (A2) and (A4) are non-linear, distribution functions (A6) 
can be only empirically found from a computer simulation of Gaussian noises 
with a given correlation function (given spectral density). 
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Figure Captions 



Fig. |1|. The stationary distribution function for binary neutron stars in the 
Galaxy dN/df /de calculated in Ignatiev et al. (2001) for Maxwellian kick 
velocities with the mean value 200 km/s. The normalization to the galactic 
binary NS coalescence rate 10~^ per year. Shwn is a fragment for systems 
with e > 0.3 which coalesce over the time interval longer than 1 year (the 
upper boundary on frequencies and eccentricities) and produce a GW am- 
plitude in the periastron h^ax estimated from Eqn (|). The signal-to-noise 
ratio is 5v^ for one sequence of pulses (|). 

Fig. ^. Upper panel: the polar angle of an elliptical orbit inside which 
99% (101 < 099) and 90% (|0| < 09o) of energy is emitted in the periastron 
passage, as a function of e. Middle panel: the pulse duration r99 and t^q 
in units of the orbital period P. Bottom panel: the variation of the mean 
amplitude h{(f))/hmax within the boundary points of the orbit 09o and ^99. 

Fig.. H. The distribution of binary NSs with GW amplitudes at periastron 
exceeding the LISA noise by 5a/5 and 5 times (lower and upper curves, 
respectively), over the pulse duration tqo (the dashed curves) and orbital 
period P (the solid curves). The stationary distribution function form Fig. 
|l] is used. The integral over the distribution for S/N = 5-\/5 (lower curves) 
is ~ 7 and for S/N = 5 (upper curves) is ~ 15. 
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Figure 2: 
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